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ABSTRACT 

Atmospheric  radio-noise  samples  were  recorded  near 
Boulder,  Colorado,  and  analyzed  to  obtain  amplitude  and  time 
statistics  of  MF  and  HF  noise.  Amplitude  and  time  statistics 
derived  from  the  HALL  model  were  compared  with  the  corresponding 
measured  statistics  for  atmospheric  noise.  The  measured  dis¬ 
tributions  of  noise  samples  previously  used  in  error-rate  tests 
were  compared  with  distributions  derived  from  the  HALL  model 
and  from  standard  atmospheric  iudio-noise  distributions.  Per¬ 
formance  curves  were  computed  for  a  CPSK  system  and  for  an 
optimum  receiver  operating  under  the  same  conditions.  Atmos¬ 
pheric  radio-noise  predictions  for  20  CONUS  locations  were 
furnished. 
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SAMSO  PHASE  C  -  FINAL  REPORT 
NOISE  DATA  AND  ANALYSIS 


A.  D.  Spaulding,  R.  T.  Disney,  and  L.  R.  Espeland 

Atmospheric  radio -noise  samples  were  recorded  near 
Boulder,  Colorado,  and  analyzed  to  obtain  amplitude  and  time 
statistics  of  MF  and  HF  noise.  Amplitude  and  time  statistics 
derived  from  the  Hall  model  were  compared  with  the  corres¬ 
ponding  measured  statistics  for  atmospheric  ncise.  The  meas¬ 
ured  distributions  of  noise  samples  used  in  previous  error- 
rate  tests  were  compared  with  distributions  derived  from  the 
Hall  model  and  from  standard  atmospheric  radio -noise  distri¬ 
butions.  Performance  curves  were  computed  for  a  CPSK 
system  and  for  an  optimum  receiver  operating  under  the  same 
conditions.  Atmospheric  radio-noise  predictions  for  20  CONUS 
locations  were  furnished. 

1.  INTRODUCTION 

The  work  covered  by  this  report  consists  of  that  performed  in 
compliance  with  task  2.10e  and  f  of  Addendum  No.  2  (Phase  C)  to 
Contract  F04701-68-F-0072,  Communications  Field  Test  FY-68, 

Project  672A  in  support  of  HRSD  Program  125B.  Task  2.10e  was  to 
record  atmospheric  radio  noise  at  suitable  times  near  Boulder,  Colorado, 
at  frequencies  of  495  kHz  and  2.  |>  MHz  and  analyze  the  recorded  noise 
samples  to  obtain  the  amplitude  probability  distributions  (APD)  and  pulse 
duration  distributions  (PDD),  pulse  spacing  (interval  between  pulses) 
distributions  (PSD),  and  average  crossing  rates  at  various  levels.  The 
recording  system  used  made  additional  recordings  possible  at  a  negli¬ 
gible  increase  in  cost.  Because  three  frequencies  at  a  time  can  be 
recorded  on  a  single  seven-track  tape,  recordings  were  also  made  at 
250  kHz  and  5  MHz.  The  requirements  of  Task  2.10f  were  to  (1)  deter¬ 
mine  the  applicability  of  the  Hall  model  for  impulsive  phenomena  to  HF 
atmospheric  radio  noise,  (2)  analytically  develop  an  optimum  (minimum 
probability  of  bit  error)  receiver  for  HF  atmospheric  radio  noise. 


(3)  Compute  the  error  rate  for  the  optimum  receiver  and  for  a  matched 
filter  linear  receiver  for  NCFSK  and  DCPSK  systems  operating  in  the 
same  noise  environment,  and  (4)  compute  the  error  rate  for  a  CPSK, 

8  bps,  4  kHz  bandwidth  system  and  the  optimum  receiver  performance 
for  the  atmospheric  radio  noise  samples  used  in  earlier  Minuteman  tests, 
comparing  the  results  with  the  earlier  results.  Much  of  the  original 
work  on  the  definition  and  specification  of  the  mathematical  modeling  of 
the  noise  was  performed  on  a  different  project  (Disney  and  Spaulding,  1968). 
Part  of  that  project  report  dealing  with  the  Hall  model  for  impulsive 
phenomena  and  its  applicability  to  atmospheric  and  man-made  noise  and 
the  derivation  and  analysis  of  performance  of  the  optimum  receiving 
system  is  appended  to  this  report  for  clarity  and  background  information. 

In  addition  to  these  requirements,  additional  analysis  of  the  atmos¬ 
pheric  radio  noise  tapes  recorded  under  Task  2.10e  provided  other  sta¬ 
tistics  of  the  noise.  The  spectrum  and  autocovariance  were  computed 
for  several  200  ms  random  segments  of  the  recordings. 

Predictions  of  atmospheric  radio  noise  were  provided  in  terms  of 
time -block  medians,  decile  deviations  from  the  medians,  and  standard 
deviations  of  these  values  within  the  time  block  for  20  possible  installation 
locations  in  the  continental  United  States. 

The  results  obtained  and  the  methods  used  in  the  above  tasks  are 
presented  in  the  following  sections. 

2.  ATMOSPHERIC  RADIO -NOISE  RECORDINGS  AND  ANALYSIS 

Estimates  of  the  amplitude  statistics  of  atmospheric  radio  noise 
are  given  in  a  report  by  CCIR  (1964),  which  does  not  contain  a  description 
of  the  time  statistics  at  various  amplitudes;  nor  are  adequate  estimates 
of  these  statistics  of  the  noise  envelope  available  elsewhere. 

Time  and  amplitude  statistics  of  atmospheric  radio  noise  in  the 
MF  and  lower  HF  portions  of  the  spectrum  are  needed  to  verify  the 
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accuracy  of  the  mathematical  model  of  the  noise.  This  information  was 
needed  also  for  theoretical  prediction  of  system  operation  when  such 
noise  is  the  predominant  type.  Because  of  time  limitations,  noise 
recordings  from  which  the  statistics  could  be  obtained  had  to  be  made 
during  a  season  of  the  year  when  daytime  atmospheric  noise  was  low  at 
the  frequencies  of  interest.  The  recordings,  therefore,  had  to  be  made 
at  night  when  atmospheric  radio  noise  was  higher  than  other  types  of 
noise.  To  obtain  some  variation  of  the  statistics,  they  were  made  on 
four  separate  nights.  One  tape  with  three  simultaneous  frequencies 
was  recorded  each  night.  The  500  kHz  and  2.  5  MHz  channels  were 
recorded  on  each  of  the  four  nights,  the  250  kHz  channel  was  recorded 
on  three  nights,  and  the  5.  0  MHz  channel  on  one  night.  The  recordings 
were  all  made  at  a  site  remote  from  power  and  telephone  lines  and  from 
highways,  and  all  were  recorded  shortly  after  midnight.  The  tape 
recording  system  described  in  the  Phase  B  final  report  (Disney  and 
Spaulding,  1968)  was  used.  The  noise  envelope  was  obtained  from  the 
analog  tapes  and  digitized  for  computer  analysis.  The  A  to  D  conversion 
was  done  at  a  10  kHz  sampling  rate.  Seven  levels  and  11  time  intervals 
were  used  in  the  computer  analysis.  The  time  intervals  throughout  were 
0.2,  0.6,  2,  6,  20,  60,  200,  600,  2000,  6000,  and  20000  ms.  The  levels 
were  limited  by  the  dynamic  range  of  the  A  to  D  converter.  A  60  dB 
dynamic  range  was  available,  and  only  the  top  60  dB  of  the  90  dB  recorded 
on  the  analog  tape  was  used.  Six  levels  10  dB  apart  (50  dB)  and  one  step 
of  8  dB  at  the  top  were  used.  The  table  1  calibration  sheet  gives  the  levels 
in  terms  of  d3  above  kT^b  for  each  of  the  following  12  computer  print¬ 
outs.  Tables  2  through  13  are  the  computer  printouts,  all  in  the  same 
format.  Tables  2,  3,  4,  and  5  are  the  results  of  the  analysis  of  the 
2.  5  MHz  recordings  for  March  29  and  April  3,  4,  and  5,  respectively; 
tables  6,  7,  8,  and  9  are  from  the  500  kHz  channel  for  the  same  respec¬ 
tive  days;  tables  10,  11,  and  12  are  from  the  250  kHz  recordings  on 
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March  29,  April  3,  and  April  5,  and  table  13  shows  the  results  of  the 
analysis  of  the  5  MHz  recording  made  on  April  4. 

The  APD's  are  given  at  the  top  of  the  page.  The  top  line  is  the 
percent  of  the  time  that  level  1  is  exceeded,  the  second  line  is  the 
percent  level  2  is  exceeded,  etc.  ,  through  level  7  on  line  seven. 

The  cumulative  PSD's  are  given  in  the  next  group  of  figures. 

The  first  column  is  the  distribution  for  level  1,  the  second  column  for 
level  2,  etc., through  the  seventh  column  for  level  7.  The  first  line  is 
the  percentage  of  the  intervals  betweer  pulses  that  are  0.  2  ms  or  more 
in  duration.  The  ether  lines  in  order  are  the  percentage  of  intervals 
that  equal  or  exceed  0,  6,  2,  6,  20,  60,  200,  600,  2000,  6000,  and 
20000  ms. 

The  next  group  of  numbers  gives  the  percentage  of  pulse  durations 
that  equal  or  exceed  the  various  lengths  of  time  in  the  same  order  as  the 
intervals  are  given.  That  is,  the  first  column  gives  the  percentage  of 
pulse  durations  that  equal  or  exceed  0.2,  0.6,  2,  6,  etc.,  ms  at  level  1. 

The  bottom  line  on  the  printout  is  the  average  crossing  rate  at 
each  of  the  levels.  The  number  of  positive  crossings  per  second  at 
level  1  through  7  are  given  from  left  to  right. 

Computer  printouts  of  the  same  information  for  man-made  radio 
noise  have  been  furnished  in  an  earlier  report  (ESSA,  1968). 

Sample  plots  of  the  data  given  in  the  tables  are  shown  on  figures  1 
through  12.  Figures  1  through  4  are  the  plots  of  the  data  given  in  table  2 
for  atmospheric  radio  noise  at  2.  5  MHz.  Figure  1  is  a  plot  of  the  ampli¬ 
tude  probability  distribution,  figure  2  presents  the  curves  of  the  PSD's, 
figure  3  shows  the  PDD's,  and  figure  4  is  a  plot  of  average  crossing 
rates.  These  values  were  recorded  on  March  29.  The  same  sets  of 
values  for  500  kHz  recorded  on  April  3  and  tabulated  In  table  7  are  shown 
on  figures  5  through  8,  and  the  data  for  250  kHz  recorded  on  April  5 
(table  12)  are  presented  in  figures  9  through  12. 
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Each  tape  was  further  analyzed  to  obtain  the  autocovariance 
function.  From  each  tape,  three  intervals  of  2048  samples  each  were 
chosen  at  random.  The  2048  data  points  correspond  to  a  time  interval 
of  essentially  0.  2  s  (10000  kHz  sampling  rate).  For  each  set  of  data, 
the  d-c  term  was  removed,  and  the  conditioned  data  were  fast-Fourier  - 
transformed  to  obtain  ar.  estimate  of  the  spectrum.  Each  data  point  is 
considered  to  be  complex,  and  the  resulting  spectrum  is  used  to  compute 
the  power  spectrum.  The  power  spectrum  (properly  scaled)  is  inverse 
fast -Fourier-transformed  to  obtain  the  autocovariance.  This  procedure 
is,  by  far,  faster  for  obtaining  the  autocovariance  than  the  standard, 
average -lagged  products  method. 

Figures  13  through  18  show  two  such  sets  of  calculations  for 
500  kHz  recorded  on  April  3.  Figures  19  through  21  show  a  set  for  500  kHz 
recorded  on  April  4,  and  figures  22  through  24  show  a  set  for  2.  5  MHz 
recorded  on  April  4.  Each  set  of  three  figures  first  shows  the  input 
(digitized)  data,  then  the  amplitude  spectrum,  and  third,  the  resulting 
autocovariance.  All  levels  given  are  in  millivolts  and  can  be  converted 
to  absolute  levels  using  table  1.  While  all  possible  (1024)  Fourier  coeffi¬ 
cients  are  computed  in  the  above,  only  the  first  50  are  given  in  the 
figures.  Each  point  shown  on  the  spectrum  corresponds  to  harmonics 
spaced  about  five  cycles  apart,  and  each  harmonic  has  its  corresponding 
negative  frequency  component.  That  is,  the  result  is  a  two-sided  spec¬ 
trum,  of  which  only  the  first  50  harmonics  of  the  positive  half  are  shown. 

Because  the  desired  result  is  the  autocovariance,  the  above  pro¬ 
cedure  does  not  give  a  particularly  good  estimate  of  the  spectrum  itself, 
since  no  conditioning  (windowing)  of  the  data  was  used.  The  spectral 
estimates  obtained  above  can  be  converted  to  a  reasonable  estimate  of 
the  true  spectrum  by  standard  conditioning  (Hamming,  for  example)  of 
these  estimates. 
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Each  of  the  above  samples  is  taken  from  within  a  burst  of  noise. 

Two  low-level  and  two  relatively  high-level  bursts  are  shown.  Figure  22 
also  shows  the  start  of  a  short  burst  at  2.  5  MHz.  Between  such  bursts 
are  periods  of  much  lower -level  background  noise. 

3.  COMPARISON  OF  MEASUREMENTS  WITH  THE  HALL  MODEL 
In  this  section  we  will  compare  the  measured  statistics  (APD, 
average  crossing  ’  ates,  PDD's,,  and  PSD's)  obtained  for  atmospheric 
radio  noise  with  the  corresponding  statistics  calculated  from  the  Hall 
model.  The  model  has  been  specified  in  a  previous  report  to  SAMSO 
(Disney  and  Spaulding,  1968),  and  the  portion  of  that  report  concerning 
the  Hall  model  and  optimum  receiver  performance  is  included  here  for 
completeness  as  appendix  A.  While  measured  APD's  and  average 
crossing  rates  have  been  compared  before  for  LF  atmospherics  (Hall,  1966), 
the  above  four  statistics  have  never  been  obtained  for  the  same  sample  of 
noise  and,  of  course,  never  compared  with  any  model. 

The  computations  of  the  envelope  and  phase  probability  distribu¬ 
tions  are  given  in  appendix  B.  Note  that  the  phase  is  uniformly  distribu¬ 
ted  and  independent  of  the  envelope,  as  one  would  intuitively  expect. 

The  APD  used  in  the  comparisons  is  given  by  (B-15) . 

The  computations  of  the  average  rate  of  envelope -level  crossings 
is  given  in  appendix  C.  Some  of  the  calculations  in  appendices  B  and  C 
were  performed  earlier  under  a  different  contract  (USAFSS  Memorandum 
of  Agreement  TT-67-334).  As  shown,  the  average  crossing  rate  is 
dependent  on  the  spectral  moments  of  the  measurement  receiver  bandpass 
characteristic  and  the  autocovariance  function  of  the  noise  process.  For 
typical  receiver  characteristics,  these  parameters  can  be  included  in 

the  rms  bandwidth,  P  ,  as  defined  in  (C-13).  Our  measurements  were 

c 

made  in  a  bandwidth  of  4  kHz,  which  corresponds  for  typical  bandpass 
shapes  to  an  rms  bandwidth,  Bc»  somewhere  between  800  ar.d  1500  Hz. 


We  will  use  1000  Hz  for  in  the  following  comparisons  so  that  the 
average  number  of  crossings  by  the  envelope  of  the  level  V  ,  N(Vq)  *s 
given  by  (C-14).  Since  we  obtained  only  the  positive  going  crossings 
from  the  tape  recordings  of  the  noise,  we  require  N(Vo)/2. 

The  determination  of  the  distribution  of  envelope -level  crossings 

(PDD's  and  PSD's)  for  a  random  process  is,  in  general,  a  classical 

unsolved  problem --even  for  Gaussian  processes.  Since  we  are  dealing 

with  a  bandlimited  envelope,  a  useful  representation  of  the  envelope  can 

be  given  in  terms  of  2TqB  samples,  along  with  a  set  of  appropriate 

orthonormal  basis  functions  (see,  for  example,  Hannock  and  Wintz,  1966, 

App.  II),  where  2B  is  the  bandwidth  and  Tq  is  the  time  interval  of 

interest.  That  is,  the  envelope  for  the  time  interval  is  given  by 

2  T  B 
o 

V(t)  =  £  ,  (1) 

i  =  1 

where  V(t.)  denotes  the  value  of  the  envelope  at  time  t^  and  4>^(t)  denotes 
the  basis  function.  Of  course,  each  V(t^)  is  a  random  variable.  Since 
the  expected  shape  of  the  received  atmospherics  is  determined  by  the 
bandpass  of  the  receiver,  $.(t)  should  be  related  to  this  bandpass  charac¬ 
teristic.  Although  generally  the  best  approximation  comes  from  using 
prolate  spheroidal -wave  functions  for  <j>^(t),  we  will  follow  Hall's  (1966) 
example  and  use  pulse -basis  functions.  The  ith -basis  functions  will  be 
a  rectangle  pulse  of  width  At  from  t^  -  At/2  to  t.  +  At/ 2  and  a  height 
1  / \J At  ,  At  denoting  the  time  interval  between  samples.  Letting  the 
number  of  samples  2TqB  be  denoted  by  N,  and}from  (C-2),  we  have  the 
envelope  given  by 

N 

V(t)  -  ^  |a(t.){  E(t.)  4>.(t)  <sT5t  .  (2) 

i  =  0 


-7- 


where,  as  in  appendix  C,  E(t)  is  the  envelope  of  the  narrow -band 
Gaussian  process,  n(t). 

We  can  now  proceed  to  compute  the  distribution  of  the  pulse 
spacings  at  the  envelope  level  V  .  Let  Prob  [T  >  Tq]  denote  the  proba¬ 
bility  that  the  time  between  a  down  crossing  of  level  Vq  to  the  next  up 
crossing,  T,  is  greater  than  T  .  Then 


Prob  [T  >  T  ]  =  Prob  frio  up  crossing  of  in  [  t,  t  +  T^]  j 
a  down  crossing  at  time  t}  . 


(3) 


The  vertical  line  in  (3)  denotes  a  conditional  probability,  i.e. ,  (3)  is  the 
probability  the  first  statement  is  true,  given  that  the  second  is  true, 
and  is,  therefore,  given  by 


Prob  [  T  >  T  ]  = 
L  oJ 


Prob  {no  up  crossing  of  Vq  in  [t»  t  +  TQ] ,  and  £ 

Prob  {down  crossing  at  t 


ana  a  aowncrossing 


In  the  range  [t»  t  +  T^J  we  have  N  samples,  where  N  =  2TqB  (i.  e. , 
smallest  integer  >  2T  B).  Thus  (4)  becomes,  in  terms  of  these  samples. 


Prob[T>TQ]  = 


Prob  [  V(tQ)  >  Vq,  Vltj)  <  Vo,  V(t2)  <  V{ 


VV  5  V 


Prob  [V(t  )  >V  ,  V(t,)  <  V  ] 
o  o  1  —  o 


If  atmospheric  noise  were  composed  of  independently  occurring 
events  (e.  g. ,  events  occurring  according  to  a  Poisson  distribution),  the 
pulse  spacings  would  be  essentially  exponentially  distributed,  especially 
at  the  higher  envelope  levels.  Some  slight  deviation  might  be  expected 
due  to  the  receiver  response  characteristic.  The  measurements  given 
indicate  that  we  are  not  dealing  with  independently  occurring  events, 
since  an  exponential  distribution  plots  as  a  straight  line  of  slope  -1  on 
the  coordinates  used  for  plotting  the  PSD's  and  FDD's  (Rayleigh  paper). 
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We  see  from  the  measured  statistics  (figs.  2,  6,  and  10) 
that  for  time  intervals  on  the  order  of  20  ms  or  less  the  PSD's 
are  reasonably  close  to  being  exponentially  distributed.  These 
time  intervals  correspond  to  events  within  a  burst  of  noise,  and 
we  have  seen  that  the  autocovariance  function  generally  shows  only 
low  correlation  for  these  short  time  periods  (figs.  15  and  24). 
However,  a  great  deal  of  correlation  exists  at  longer  time  periods 
(time  between  bursts,  for  example),  as  reflected  by  the  much 
greater  steepness  of  the  PSD's  at  the  longer  time  intervals  and 
higher  envelope  levels. 

We  will  calculate  (5)  and  the  corresponding  expression  for 
PDD's  for  two  cases.  In  the  first  case  we  will  assume  independent 
samples,  and  in  the  second  we  will  introduce  some  correlation 
between  samples. 

If  the  samples  are  independent,  then  (5)  becomes 


Prob  [  T  >  T  ]  = 


Prob  [  V(t)  >  V  ]  (Prob  [  V(t)  <  V  ]f 
_  o _ —  o 

Prob  [  VCt)  >  Vq]  Prob  [  V(t)  <  V  ] 


so  that 


N  -  1 

Prob[  T  >  T  J  =  {Prob  [V(t)  <  V  1} 
o  —  o 


and,frorr.  the  expression  for  the  envelope  distribution  (B-15), 


(6) 


(7) 


Prob[  T  > 


m 


G/vo) 

♦  (v/vj2_ 


N 


1 


(8) 


The  result  (8)  corresponds  to  one  previously  obtained  by  Hall  (1966). 
Similarly,  if  we  let  Uq  denote  the  time  interval  of  interest,  then  for  the 
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pulse  durations,  i.e.  ,  the  time  between  an  up-going  crossing  and  the 
next  down-going  crossing  of  level  V  ,  we  have 


Prob  [  U  >  U  I  = 


Prob  [  V(t  )  <  V  ,  V(t  )  >  V  ,  V(t  )  >  V . V(t  )  >  V  ] 

o  o  l—o  2—o  N  —  o 


Prob  [  V(t  )  <  V  ,  V(t.)  >  V  ] 
o  o  1  _  oJ 


(9) 


which  reduces,  for  independent  samples,  to 

N  -  1 

Prob  [  U  >  U  ]  =  {Prob  [  V(t)  >  V  ])  ,  (10) 

o  **■  o 

so  that 

N  -  1 

f  (Y/Vo)m  1 

Prob[U>Uo]=  | (11) 

[l  +  (Y/Vo)2] 

We  next  introduce  some  correlation  into  the  noise  process. 

This  means  that  a  covariance  function  must  be  specified  for  the  process  a(t)^ 

(The  independence  assumption  above  corresponds  to  specifying  a  delta 

function  covariance. )  We  will  let  a(t)  remain  constant  over  the  time- 

interval  of  interest  (T  or  U  ),  so  that  a(t)  is  now  a  random  variable 

o  o 

rather  than  a  random  process. 

Let  P_t(V  )  denote  the  numerator  of  (5), 

N  o 

Prob  [  V(t  )  >  V  ,  V(t_)  <  V  ,  V(t  )  <  V . V(t  )  <  Vl  . 

o  o  i  —  o  2—o  In  —  o 

W  e  now  have 


PN(Vo>  -  J  dV0  J  dVl  • ' ' '  J  dVN  Pv  (  ’  <l2> 

Vo  o 

o 
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where  V  denotes  V(t.),  and  V  denotes  the  vector  of  samples 
1  i  — 

{Vq,  Vj . V^}.  To  maintain  notation  consistent  with  previous  work, 

we  now  have  an  ambiguity.  Although  Vq  is  the  particular  envelope  level 

of  interest,  in  (12)  V  is  also  used  for  the  sample  V(t  ).  The  context 

o  o 

is  such,  however,  that  no  confusion  will  arise. 

Since  our  envelope  V(t)  over  the  time  interval  Tq  is  now  given 
by  ( see  app.  C) 


V(t)  =  i  a  I  E(t)  , 


(13) 


where  E(t)  is  the  envelope  of  the  narrow-band  Gaussian  process  n(t), 


Py  (V)  =  |  -fVr  P|  a  |  (x)  »E  (=)  • 


(14) 


where  E  denotes  the  vector  {E. }  of  N  samples  corresponding  to  the 
vector  {V.}.  Now  p  (x)  is  given  by  (C-7)  and,  since  E  is  the  envelope 

1  lal 

of  a  Gaussian  process,  each  E.  is  Rayleigh  distributed  and  independent 


so  that  (see  app.  A) 


N 


p  (E)  =  n 
±  -  i  =  0 


(15) 


The  integral  (14)  can  now  be  evaluated,  giving 

/  2N  +  m  +  2  \ 

■+i  m±ZZZZZ _ 

r(r) 


PV(V)  =  2 


N 

n  v. 

i=  o  1 

2N  +  m  +  t  * 


N 

I 

.i=  1 


2  2 
+  y 


(16) 
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Substituting  (16)  into  (12)  and  evaluating  the  N  +  1  dimensional  integral, 
we  have 


pM(V  )  =  \ 
N  o 


m 


N 

I 


,  „k(*0 

(-D  \k  / 

r  2 

?1m/?. 

[(k  +  1)  vo  + 

J 

(17) 


Likewise,  for  pulse  durations  (time  interval  U  ),  P  ( y  )  is  now  given  by 

o  No 


W  =  J  dVo  fdVl-"  |dVNPV^» 


(18) 


which  gives 


PN(Vo>  * 


r  2 

zim/2  r  ?  2-i 

(N+l)  Vq  +  y 

(N  +  2)  V  +  y 

L  0 

m 


72  * 


Finally  then,  for  this  case  of  dependent  samples,  (5)  gives 


Prob  [  T  >  T  ]  =  — - 


k  +  1  +  (y/vq) 


~\m/Z 


(-D 


Z  f  7  "|m/2 

=  0  k  +  1  +  (y/V  ) 

L  O 


i.nd  (9 )  gives 


Prob  [U  >  U  ]  = 
o 


I 


(N  +  1)(Vo/y)2  + 1] 


m/2 


[(N  + 


(19) 


(20) 


2)(Vo/yr  +  1 


r 


72 


(21) 


r  2  I  nn/2 

r  2  i 

[(voA)  +ij 

[2(Vo/v>  +Ij 

m 


72 
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Figures  25  through  34  show  the  comparison  of  the  measured 

envelope  distributions  and  average  crossing  rates  with  those  calculated 

(B-15)  and  (C-14)  for  the  Hall  model  for  five  of  the  atmospheric -noise 

tapes  analyzed,  two  at  500  kHz,  two  at  2. 5  MHz, and  one  at  250  kHz.  The 

levels  are  given  relative  to  the  parameter  y.  since  the  theoretical  results 

are  in  terms  of  V  /y.  They  can  be  converted  easily  to  absolute  values 

(dB  >  kT  b),  as  explained  earlier.  The  comparison  shows  good  agree  - 
o 

ment.  and  for  the  average  crossing  rates  could  be  made  even  better  by 
specifying  the  proper  spectral  moments  (app,  C),  r  ither  than  the  "loose" 
parameter  B^. 

Figures  35  and  2  3  show  the  comparison  of  the  measured  PSD's 
an  1  FDD's  with  those  derived  from  the  Hall  model,  considering  independent 
samples  ((8)  and  (11)).  The  comparison  is  shown  only  for  the  250  kHz  case. 
We  see  that  the  derived  distributions  are  essentially  exponential  espe¬ 
cially  at  the  low  probabilities,  and  do  not  match  the  measured  distribu¬ 
tions  for  the  larger  time  intervals. 

Figures  37  and  38  show  fhe  comparison  for  considering  dependent 
samples  ((20)  and  (21)).  which  indicates  too  much  dependence  at  the 
larger  time  intervals  as  characterized  by  the  derived  distributions 
becoming  too  steep  and  grouping  together  more  rapidly  than  the  meas¬ 
ured  distributions  at  the  higher  envelope  levels. 

While  the  amplitude  statistics  and  the  average  time  statistics 
(average  number  of  crossings  per  second)  match  quite  well,  the  FDD's 
anl  PSD's  do  not  matcn  as  well.  The  two  cases  given  (independent 
samples  and  some  degree  of  correlation)  bracket  the  real  situation.  It 
is  apparent  that  a  proper  covariance  function  can  be  specified  for  a(t), 
so  that  these  time  statistics  also  can  be  closely  matched. 

Note  that  (2)  is  quite  difficult  to  evaluate  if  N  is  at  all  large 
because  of  the  alternating  series,  coupled  with  the  binomial  coefficients. 
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.  For  example,  in  our  4  kHz  bandwidth,  a  time  interval  of  2000  ms 
corresponds  tc  an  N  of  8000  Techniques  have  been  developed  to  evaluate 
(20)  and  computer  programs  written  to  implement  them. 

Consider  also  the  special  case  0  =  3(m  =  2),  where  (20)  can  be 
evaluated  easily,  since 

y  (-» kM 

k=o  [k*1+<y/V0>2] 

and  the  integral  can  be  readily  evaluated  numerically. 

4.  CALCULATION  OF  PERFORMANCE  OF  OPTIMUM  RECEIVER 
FOP  CURRENT  MINUTEMAN  SYSTEM 

The  current  Minuteman  receiver  uses  limiting  of  the  noise  in  a 

4  kHz  bandwidth,  followed  by  an  effective  bandwidth  reduction  to  8  Hz. 

The  receiver  is  a  standard  DCPSK  receiver  operating  at  8  bps.  We  will 
apply  the  results  of  appendix  A  to  the  above  situation  to  obtain  the  per¬ 
formance  of  the  optimum  receiver. 

Many  tests  of  the  effect  of  limiting  and  hole  punching  on  the  per¬ 
formance  of  the  above  standard  receiver  have  been  conducted  (Sisco, 

1964,  Sylvania,  1963),  and  the  object  here  is  to  determine  how  much 
better  the  optimum  receiver  would  perform.  Since  atmospheric -noise 
tapes  were  used  in  the  tests  cited,  we  must  first  specify  the  parameters 
of  the  Hall  model  that  match  the  statistics  of  the  noise  used.  The  noise 
had  essentially  the  values  of  V^(app.  A)  of  26.  0  dB,  12.. 0  dB,  and  7.  0  dB. 
Figure  39  i8  a  comparison  of  the  measured  distribution  and  the 

distribution  from  the  Hall  model  for  m  =  1.  0  (V ,  =  26.  0  dB)  for  tracks 

d 

5  and  7,  reel  1,  of  the  Canal  Zone  tape.  Figure  40  shows  the  standard 
atmospheric -noi 8 e  distribution  for  =  7.  0  dB  ,  with  the  distribution 


f 


(\/V 


o) 


(1  -  z)N  dz 


(22) 
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from  the  Hall  model  for  m  =  2.  5.  Figure  41  is  a  comparison  of  the 
measured  distribution  from  tracks  1  and  3,  reel  2,  of  a  Canal  Zone 
tape,  along  with  distributions  for  =  12.  0  dB  and  m  =  1.  ,7.  Figure  42 
shows  the  measured  distribution  from  tracks  1  and  3,  reel  1,  Canal  Zone 
tape  along  with  the  distribution  for  m  =  1.  0.  The  value  =  26.  0  dB 
would  correspond  here  also.  The  frequency  throughout  is  450  kHz,  and 
since  the  bandwidth  is  4  kHz  we  require  the  performance  of  the  optimum 
receiving  system  for  m  =  1.  0,  1.7,  and  2.  5. 

The  tests  reported  by  Sisco  (1964)  and  Sylvania  (1963)  were  generally 
conducted  with  the  receiver  operating  in  a  CPSK  mode  rather  than  DCPSK. 
We  will  therefore  compute  performance  for  CPSK  signale. 

First,  for  comparison,  we  wall  give  the  performance  of  the 
standard  CPSK  receiver  operating  at  8  bps.  The  4  kHz  noise  distribu¬ 
tions  of  V  ,'s  of  26.  0,  12.  0,  and  7.  0  dE  reduce  in  an  8  Hz  bandwidth  to 
d 

distributions  with  V^'s  of  3.  7,  1.  83,  and  1.  55  dB  (Spaulding  et  al.  ,  1962). 
The  performance  for  these  three  values  of  is  shown  in  figure  43  and 
were  computed  by  Haiton  and  Spaulding  (1966)  and  Spaulding.  Also 
shown  is  the  performance  curve  for  the  standard  (linear,  matched  filter 
CPSK  receiver  in  Gaussian  noise. 

Now  it  remains  to  compute  the  performance  of  the  optimum 
receiving  system  for  the  corresponding  situations.  As  shown  in  (A-35), 
performance  (probability  of  error,  p^)  is  given  by 


2TB 


p  =  Prob 
e 


IS 


In  z.  > 


(23) 


where  2B  is  the  bandwidth,  T  is  the  time  duration  of  the  signal,  and  the 
random  variables  z^  are  functions  of  the  complex  envelopes  of  the  noise 
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and  signal,  as  given  in  (A-36),  In  our  situation  then,  2TB  is  500..  In 
appendix  A,  performance  was  calculated  for  2TB  equal  10  and  20,  and 
was  given  by  the  appropriate  Chernoff  upper  bound  (A -46).  Here  the  high 
value  of  2TB  results  in  this  upper  bound  not,  in  itself,  being  a  very  good 
estimate,  which  must  be  modified  according  to  (A -41).  That  is,  the 
factor  i/sl  2n  s2  •  •  must  be  calculated.  While  this  factor  was  of 

little  significance  in  the  previous  calculation,  it  now  becomes  highly 
significant,  as  can  be  seen  from  (A-52),  which  gives  ‘p'(s)  for  the  special 
case  m  =  2.  Numerical  techniques  were  developed  to  evaluate  p‘ (s)  for 
the  values  of  m  of  interest  here.  Even  so,  since  we  are  dealing  with  low 
signal -to -noise  ratios,  the  resulting  estimate  (A-41)  is  valid  only  for 
quite  low  pe> 

To  complete  the  performance  estimates,  we  now  will  obtain  a 
lower  bound  that  is  quite  good  for  low  signal -to -noise  ratios  and  high 

Pa¬ 
using  the  inequality, as  suggested  by  Hall  (1966), 


we  have 


In  z.  >  1 

l  — 


_1_ 

z. 

l 


9 


Pe  — 


> 


(24) 


(25) 


Since  2TB  is  500,  we  now  apply  the  central  limit  theorem  to  (25),  that 
is,  we  require  the  mean  and  variance  of  the  random  variables  1  -  l/z.. 


I 

I 

I 


F rom  appendix  A,  we  have 


1  - 


_i_ 

z. 

1 


(26) 


where  h.  is  the  ith  sample  of  the  complex  noise  envelope,  P.  is  the  ith 

1  2.  2  i 

sample  of  the  complex  signal  envelope,  and  %  is  mcr^  A t/2(r  (app.  A). 

With  ti  =  n  +  j  n.  and  p.  =  p.  +  j  p.  (that  is,  rj.  and  n.  are  the 
i  icis  l  ic  is  ic  is 

real  and  imaginary  parts  of  the  complex  random  variable  ^  ),  we  have 


00  00 


E[1-r]= 


00  -00 


^|2+  Mi'^il21 


\\  + 


2  .  ^ 


>  Ai.  .  h.  \ 

V’  ^is  \  C  18j* 


(27) 


where  (Hall,  1966) 


/  \  m  (; 

\e,  r'is\'ic'  T,i8)"  2Tr 


m  +  2 


;>**♦*; 


and  E  denotes  the  expectation  operation  (mean  value  operation). 
Evaluating  (27),  we  obtain 


L  V  (m  +  2) 


(28) 


(29) 


so 


H)]  ■  •  - 1  ^ 


Similarly, 


j’j  _  j_j 2  =  4m  I 


(m  +  2)(m  +  4) 


(m  +4)  £ 


(30) 


(31) 
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X 


The  variance  of  the  variable  1  -  l/z^  is  then  given  by 

2 


n  u 

4m 

“i 

i 

h 

3 

1 

3 

** 

*1 

i 

j 

r*j 

(m  +  2)(m  +  4) 

£ 

\m  +  4  <m+2)2  / 

£ 

• 

Finally,  with 


2TB  2 

T-  I  I 


H  =  -  — 

z  m 


m  V'  'i 

+  2  Z  5 


2TB  ..  2 


2  4m 

<r  =  — - — 

z  (mv4)(m 


we  obtain  an  estimate  of  the  lower  bound  given  by  (25) 

pe  i ertc  (r1)  • 


Here 


00 

trfc (x)  =  y  exp  (  ~yZ/2)  dy  . 


By  using  both  the  upper  bound  (A-41)  and  the  lower  bound  (35),  the 
performance  curves  for  m  =  1.  0,  1.  7,  and  2.  5,  given  in  figure  43, 
were  obtained.  Note  that  the  above  calculations  are  for  signals 
with  zero  crosscorrelation  (app.  A),  so  that  the  curves  of  figure  43 
also  incorporate  the  standard  3  dB  improvement  one  expects  in 
going  to  signals  with  -1  crosscorrelation  (as  in  the  case  of  CPSK 
signals). 

While  the  results  of  the  limiter  tests  (Sisco,  1964;  Sylvania, 
1963)  are  not  included  in  figure  43,  they  generally  group  about  the 
Gaussian  characteristic,  indicating  that  even  in  this  rather  extreme 
case  of  limiting  and  bandwidth  reduction,  the  optimum  receiver 
still  performs  substantially  better. 

-18- 


i  5.  REFERENCES 

t 

Akima,  H.  (1967),  Speeding  up  of  HF  radio-transmissions  and  broadcasts 

*  (World  Weather  Watch  Study  T ..  21),  ESSA  Tech.  Memo.  IERTM-ITS^  42.. 

Beckmann,  P.  (1964),  Amplitude  probability  distribution  of  atmospheric 

*  radio  noise,  Radio  Sci.  68D,  No.  6,  723-736. 

Bello,  P.A.  (1965),  E-.-ror  probabilities  due  to  atmospheric  noise  and 
|  flat  fading  in  K...r  ionospheric  communications  systems,  IEEE 

|  Trans.  Commun .  Technol.  COM -13,  No.  3,  266-279. 

\  Bowen,  B.  A.  (1963),  Some  analytical  techniques  for  a  class  of  non- 

[  Gaussian  processes.  Research  Report  No.  63-3,  Queens  University, 


CCIR  (1964),  World  distribution  and  characteristics  of  atmospheric  radio 
noise.  Report  322,  ITU,  Geneva. 

Conda,  A.  M.  (1965),  The  effects  of  atmospheric  noise  on  the  probability 
of  error  for  an  NC  'SK  system,  IEEE  Trans.  Commun  .  Technol. 

COM -13,  No.  3,  280-283. 

Crichlow,  W.  Q. ,  A.  D.  Spa  ilding,  C.  J.  Roubique,  and  R.  T.  Disney 
(I960),  Amplitude  Probability  Distributions  for  Atmospheric 
Radio  Noise,  NBS  Monograph  2?  (U.S.  Government  Printing  Office, 
Washington,  D.  C.  ). 

Disney,  R.  T.  and  A.  D.  Spaulding  (1968),  Interference  measurements  and 
analysis.  Final  Rept.  to  Spree  and  Missile  Systems  Organization,  VI, 
prepared  by  the  Environmental  Science  Services  Administration  under 
Contract  No.  F04'r01-69  0072,  Project  67 2A. 

ESSA  (1967),  Communications  Task  -  Project  672A,  Final  Rept.  to  Space 
and  Missile  Systems  Organization,  prepared  by  the  Environmental 
Science  Services  Administration  under  Contract  AFD/0  (04 -691) -67 -3. 

ESSA  (1968),  Interference  field  tests,  Interim  Evaluation  Rept.  ,  prepared 
by  the  Environmental  Science  Services  Administration  under 
Contract  No.  F04701-68-F-0072,  Project  672A. 

Furutsu,  K. ,  and  T.  Ishida  (I960),  On  the  theory  of  amplitude  distributions 
of  impulsive  random  noise  and  its  application  to  atmospheric 
noise,  Japan  Radio  Res.  Labs.  J.  ?,  No.  32,  279-307. 


-19- 


Galejs,  J.  (1966),  Amplitude  distributions  o£  radio  noise  at  ELF  and 
VLF,  J.  Geophys.  Res.  71,  201-216. 

Hall,  H.  M.  (1966),  A  new  model  for  "impulsive"  phenomena:  application 
to  atmospheric -noise  communication  channels.  Tech.  Rept. 

No.  34l2-8  and  7050-7,  Stanford  Electronic  Laboratories, 

Stanford,  California. 

Halton,  J.  H. ,  and  A.  D.  Spaulding  (1966),  Error  rates  in  differentially 
coherent  phase  systems  in  non -Gaussian  noise,  IEEE  Trans. 
Commun-  Technol.  COM -14,  No.  5,  594-601. 

Hannock,  J.  C.  and  P.  A.  Wintz  (1966),  Signal  Detection  Theory, 
(McGraw-Hill  Book  Co.  ,  Inc.,  New  York). 

Helstron,  C.  W.  (I960),  Statistical  Theory  of  Signal  Detection  (Pergamon 
Press,  Oxford,  England). 

Ibukun,  O.  (1966),  Structural  aspects  of  atmospheric  radio  noise  in  the 
tropics,  Proc.  IEEE  54,  361-367. 

Lindenalub,  J.  C.,  and  K.  A.  Chen  (1965),  Performance  of  matched 

filter  receivers  in  non-Gaussian  noise  environments,  IEEE  Trans. 
Communi.  Technol.  13,  No.  4,  545-547. 

Linfieid,  R.  F.  (1965),  High  performance  reliable  VLF  component  of 

the  Naval  Advanced  Communication  System  -  VLF  system  design 
specification  and  evaluation  process.  Final  Report  No.  34-F, 
prepared  by  DECO  Electronics,  Inc.  forU.S.  Navy  Bureau  of  Ships. 

Mandelbrot,  B.  (1964),  Self-similar  turbulence  and  non -Wienerian 
conditional  spectra  (unpublished  IBM  research  paper). 

Middleton,  D.  (1961),  Introduction  to  Statistical  Communication  Theory 
(McGraw-Hill  Book  Co. ,  Inc. ,  New  York). 

Shchukin,  A.  N.  (1946),  On  a  method  of  combating  impulse  interference 
in  radio  reception.  Bull.  Academy  of  Sci.  USSR,  Phys.  Series  10, 
No.  1,  49-56. 

Shepelavey,  B.  (1963),  Non-Gaussian  noise  in  binary  data  phase  coherent 


-20- 


Sisco,  W.  B.  (1964),  Bit  and  message  error  rate  in  atmospheric  noise 
with  and  without  peak  limiting.  Report  No.  7070-6323 -RU-000, 
TRW  Space  Technology  Laboratories,  Redondo  Beach, California. 

Spaulding,  A.  D.  ,  C.  J.  Roubique,  and  W.  Q.  Crichlow  (1962), 

Conversion  of  the  amplitude  probability  distribution  function  for 
atmospheric  radio  noise  from  one  bandwidth  to  another,  J.  Res. 
NBS  66D,  713-720. 

Spaulding,  A.  D.  (1964),  Determination  of  error  rates  for  narrowband 
communication  of  binary -coded  messages  in  atmospheric  radio 
noise.  Proc.  IEEE  52,  No.  2,  220-221. 

Spaulding,  A.  D.  (1966),  The  characteristics  of  atmospheric  radio 
noise  and  its  effects  on  digital  communications  systems, 

IEEE  Intern.  Commun.  Conf.  ,  Paper  No.  CP  1126. 

Sylvania  (1963),  Design  review  report  for  project  124-221  (Contained 
in  report  on  hole  puncher  performance  investigation  of  radio 
subsystems  of  Minuteman  Ground  Electronic  Systems),  Sylvania 
Electronics,  Waltham,  Mass. 

Van  Trees,  H.  L.  (1968),  Detection,  Estimation,  and  Modulation  Theory, 
Part  1  (John  Wiley  and  Sons,  Inc,  ,  New  York). 


-21- 


Table  1.  Calibration  Sheet 

Absolute  levels  corresponding  to  levels  1  through  7  on 
computer  printouts.  Levels  given  to  nearest  dB. 
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Table  4.  Distributions  of  atmospheric  radio  noise.  Eggleston  Reservoir 
No.  4,  Boulder,  Colorado,  April  4,  1969,  0000-0400  hrs. 
Frequency:  2.  5  MHz.  Impulse  bandwidth:  4  kHz. 
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Table  8.  Distributions  ol  atmoenhrric  radio  nois*>.  Eggleston  Reservoir 
No.  4,  Boulder,  Colorado,  April  4,  19fV,  0000-0400  hrs. 
Frequency:  500  kHz.  Impulse  bandwidth:  4  kHz. 
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Table  12:  Distributions  of  atmospheric  radio  noise.  Eggleston  Reservoir 

No.  4,  Boulder,  Colorado,  April  5,  196.9,  0000-0400  hrs. 
Frequency:  250  kHz.  Impulse  bandwidth:  4  kHz. 
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Figure  1.  Amplitude  probability  distribution  of  atmospheric  radio  noise 
recorded  at  Eggleston  Reservoir  No.  4,  Boulder.  Colorado. 
F  =  2.  5  MHz,  Bi  =  4  kHz,  March  29,  1969,  0000-0400  hrs. 
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Figure  2.  Pulse  spacing  distributions  of  atmospheric  radio  noise  recorded  at 
Eggleston  Reservoir  No.  4,  Boulder*  Colorado. 

F  =  2.  5  MHz,  B.  =  4  kHz,  March  29,  1969,  0000-0400  hrs. 
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Figure  5.  Amplitude  probability  distribution  of  atmospheric  radio  noise 
recorded  at  Eggleston  Reservoir  No.  \  x  Boulder,  Colorado. 
F  =  500  kHz,  B.  =  4  kHz,  April  3,  1969,  0000-0400  hrs  . 
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Figure  i.  Pulse  spacing  distributions  of  atmospheric  radio  .noise  recorded 
at  Eggleston  Reservoir  No.  4,.  Boulder,  Colorado. 

F  =  500  kHz,  B.  =  4  kHz,  April  3,  1969,  0000-0400  hrs. 
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Figure  7. 


rec°rded 

i  4  kHz,  April  3,  1969,  0000-0400  hrs. 


Figure  8.  Average  crossing  rates  of  atmospheric  radio  noise  recorded 
at  Eggleston  Reservoir  No.  4,  Boulder,  Colorado. 

F  =  500  kHz,  B  =  4  kHz,  April  3,  1969,  0000-0400  hrs. 
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250  kHz,  4  kHz  Bandwidth 
O  Denotes  Measured  Points 
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Figure  9.  Amplitude  probability  distribution  of  atmospheric  radio  noise 
recorded  at  Eggleston  Reservoir  No.  4,  Boulder,  Colorado. 
F  =  250  kHz,  B  =  4  kHz,  April  5,  1969,  0000-0400  hrs. 
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Figure  10.  Pulse  spacing  distributions  of  atmospheric  radio  noise 

recorded  at  Eggleston  Reservoir  No,  4,  Boulder,  Colorado. 
F  =  250  kHz,  B£  a  4  kHz,  April  5,  1969,  0000-0400  hrs .  ' 
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Figure  11.  Pulse  duration  distributions  of  atmospheric  radio  noise  recorded 
at  Eggleston  Reservoir  No.  4,  Boulder,  Colorado 
F  =  250  kHz,  B.  =  4  kHz,  April  5,  1969,  0000-0400  hrs. 
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Figure  23.  Spectrum  of  atmospheric  radio  noise  from  recording  at 
Eggleston  Reservoir  No.  4 ,  Boulder*  Colorado. 

F  a  2.  5  MHz,  B.  a  4  kHz,  April  4,  1969* 
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Figure  25.  Comparison  of  measured  amplitude  probability  distribution  with 
that  calculated  from  the  Hall  model. 
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Figure  26.  Comparison  of  measured  average  crossing  rates  with 
those  calculated  from  the  Hall  model. 
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Figure  27.  Comparison  of  measured  amplitude  probability  distribution 
with  that  calculated  from  the  Hall  model. 
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Figure  28.  Comparison  of  measured  average  crossing  rates  with  those 
calculated  from  the  Hall  model. 
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Figure  29.  Comparison  of  measured  amplitude  probability  distribution  with 
that  calculated  from  the  Hall  model. 
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Figure  30.  Comparison  of  measured  average  crossing  rates  with  those 
calculated  from  the  Hall  model. 
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Figure  31.  Comparison  of  measured  amplitude  probability  distribution  with 
that  calculated  from  the  Hall  model. 
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Figure  33.  Comparison  of  measured  amplitude  probability  distribution 
with  that  calculated  from  the  Hall  model. 
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Figure  34.  Comparison  of  measured  average  crossing  rates  with  those 
calculated  from  the  Hall  model. 
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Figure  35.  Comparison  of  measured  pulse  spacing  distributions  with 

distributions  calculated  from  the  Hall  model  (no  correlation  assumed). 
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Figure  36.  Comparison  o£  measured  pulse  durations  with  distributions 
calculated  from  the  Hall  model  (no  correlation  assumed)’. 
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Percent  of  Pulse  Spocings  Which  Exceed  Ordinate 
Figure  37.  Comparison  of  measured  pulse  spacing  distribution  with  distri¬ 
butions  calculated  from  the  Hall  model  (some  correlation  assumed). 
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Figure  38.  Comparison  of  measured  pulse  durations  with  distributions 
calculated  from  the  Hall  model  (some  correlation  assumed). 
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Figure  40. 
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Comparison  of  atmospheric  radio  noise  amplitude  probability 
distribution  for  Vd  =  7  dB,with  distribution  calculated  from  the 
Hall  model  for  m  =  2.  5. 
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Amplitude  in  dB  Relative  to  r.m.s. 


Figure  42. 


Hercent  of  Time  Ordinate  is  Exceeded 

^ri8°n  of1m®a8u  red  amplitude  probability  distribution  and 
distribution  calculated  from  the  Hall  model  for  m  *  1.  o  (V  *  26. 
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Figure  43. 


°f1the  performaace  of  ^  optimum  receiver  with  a 
matched  filter  linear  receiver  for  a  binary  CPSK  system. 
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APPENDIX  A 


A.  1  Mathematical  Modeling  of  the  Noise  Process 
To  be  able  to  determine  the  performance  of  an  optimum  receiver 
subjected  to  man-made  radio  noise,  a  model  for  this  narrow -band 
random  noise  process  must  be  developed.  Narrow -band  noise  processes 
result  whenever  the  receiver  bandwidth  is  substantially  less  than  the 
receiver  center  frequency.  The  received  waveform  of  the  narrow -band 
process  has  an  envelope  and  a  phase.  For  determining  the  optimum 
receiver,  more  information  about  the  noise  is  required  than  generally 
can  be  obtained  from  measurements.  The  problem  is  to  develop  a  model 
for  the  noise  that  fits  all  the  available  measurements,  that  is  physically 
meaningful  when  the  nature  of  the  noise  sources  are  considered,  and 
that  Btill  is  simple  enough  that  the  required  statistics  can  be  obtained 
for  solving  signal  detection  problems. 

Many  attempts  have  been  made  to  model  narrow -band  impulsive 
noise  processes  (Furutsu  and  Ishida,  I960;  Middleton,  1961;  Beckmann, 
1964;  Galejs,  1966).  These  models  are  essentially  similar  in  that  they 
take  the  received  noise  to  be  composed  of  a  sum  of  filtered  impulses 
whose  amplitudes  and  occurrence  in  time  follow  various  probability 
distributions.  For  example,  the  amplitudes  of  the  impulse  responses 
may  be  assumed  to  be  log-normally  distributed,  while  their  times  of 
occurrence  may  be  assumed  to  be  Poisson  distributed. 

Although  the  above  forms  are  well-motivated  physically  and  cam 
be  made  to  fit  measured  first-order  statistics  (amplitude  probability 
distribution  of  the  noise  envelope,  for  example),  they  have  several  dis¬ 
advantages  as  far  as  the  signad  detection  problem  is  concerned: 

(1)  The  models  assume  independence  in  the  noise,  while  meas¬ 
urements  indicate  that  this  is  not  the  case  for  either  atmospheric  noise 
or  most  forms  of  mam -made  noise.  When  one  considers- the  correlation 
in  the  sources  of  man-made  noise  (automobile  ignition  systems  resulting 
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in  essentially  periodic  pulse  trains,  for  example),  the  assumption  of 
independence  is  obviously  wrong. 

(2)  The  resulting  probability  distributions  are  quite  complicated 
and  usually  cannot  even  be  put  in  closed  form.  For  this  reason  no 
attempt  has  ever  been  made  to  apply  these  models  to  statistical  detection 
problems  at  the  receiver. 

Various  empirical  models  have  been  developed  and  related  to 
measurements  (Crichlow  et  al. ,  I960;  Spaulding  et  al. ,  1962;  Ibukun, 

1966).  These  models  do  not  represent  the  noise  process  but  only  the 
measured  statistics  of  the  process  and,  therefore,  are  not,  in  general, 
applicable  in  determining  optimal  receivers  for  the  particular  noise  under 
consideration.  They  have  been  used  to  determine  performance  of  various 
auboptimum  linear  receivers.  The  receivers  now  in  use  are  designed 
to  be  optimum  in  Gaussian  noise,  and  the  performance  of  these  receivers 
acting  in  impulsive  noise  has  been  determined  (Akima,  1967;  Bello,  1965; 
Conda,  1965;  Halton  and  Spaulding,  1966;  Lindenlaub  and  Chen,  1965; 
Shepelavey,  1963;  Spaulding,  1964,  1966)  . 

Here  we  are  considering  digital  systems  in  which  each  signaling 
element  is  equally  probable,  and  since  the  "cost"  associated  with  making 
an  error  is  the  same  for  each  type  of  possible  error,  "optimum"  here 
simply  means  minimum  probability  of  error.  For  example,  in  a  binary 
system,  by  equal  "costs  "  we  mean  making  the  error  "decide  signal  one 
when  signal  two  waa  sent"  is  just  as  bad  as  making  the  error  "decide 
signal  two  when  signal  one  was  sent. "  Indications  as  to  performance 
improvement  by  use  of  nonlinear  elements  (hole  punchers,  limiters, 
etc. )  have  also  been  given  (ESSA,  1967;  Linfield,  1965;  Shchukin,  1946; 
Sisco,  1964;  Syl vania,  1963). 

Recently,  Hall  (1966)  applied  work  on  the  applicability  of  a  class 
of  "8  elf -similar"  random  processes  as  a  model  for  certain  intermittent 
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phenomena  to  signal  detection  problems  considering  LF  atmospheric 
noise.  The  concept  introduced  is  that  of  a  random  process  controlled  by 
one  "regime"  for  the  duration  of  observation,  while  this  regime  is  itself 
a  random  process.  This  concept  was  introduced  by  Mandelbrot  (1964)  in 
the  study  of  turbulence. 

Hall  (1966)  has  shown  that  the  model  to  be  specified  in  the  next 
section  fits  the  measured  statistics  of  LF  atmospheric  noise.  Other 
work  at  ITS  has  indicated  that  the  model  may  be  modified  to  fit  HF  atmos¬ 
pherics  and  man-made  noise. 

In  the  following  section  we  will  specify  the  model  and  possible 
modifications,  determine  the  optimum  receiver  for  some  examples  of 
LF  atmospheric  noise,  considering  frequencies  around  500  kHz,  obtain 
estimates  of  this  optimum  receiver's  performance,  and  show  how  this 
performance  compares  with  the  performance  of  typical  linear  (matched 
filter  or  correlation)  receivers. 

As  expected,  the  optimum  receiver  is  nonlinear  and  its  perfor¬ 
mance  therefore  depends  on  more  than  just  the  signal -to -noise  ratio. 

We  will  see  that,  in  general,  performance  will  depend  on  the  actual 
noise  power,  the  actual  signal  energy,  the  signal  shape,  bandwidth,  and 
time  duration  of  the  signals.  We  will  analyze  our  optimum  receiver's 
performance  for  the  following  situations: 

(a)  1000Hz  bandwidth,  binary  system,  signaling  rate  of  100  bps, 
with  the  highest  atmospheric  noise  levels,  as  predicted  by 
CCIR  (1964)  for  the  northwest  U.S. 

(b)  Same  as  (a),  but  with  a  signaling  rate  of  50  bps. 

For  comparison  with  current  receivers,  case  (a)  would  correspond  to  a 
receiver  that  tries  to  fight  the  noise  by  limiting  in  a  1000  Hz  bandwidth 
and  then  uses  approximately  a  10-to-l  bandwidth  reduction  to  limit  the 
noise. 
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A.  2  Specification  of  the  Hall  Model 
The  model  proposed  for  received  impulsive  noise  is  one  that  takes 
the  received  noise  to  be  a  narrow -band  Gaussian  process  multiplied  by  a 
weighting  factor  that  varies  with  time.  We  will  consider  the  narrow -band 
noise  process  y(t),  to  be  of  the  form 

y<t)  =  a(t)  .  n(t)  ,  (A-l) 


where  n(t)  is  a  zero -mean  narrow -band  Gaussian  process  centered  on  the 

frequency  of  interest,  uo  ,  with  covariance  function,  R  (i),  and  a(t),  the 

o  n 

"regime"  process,  is  stationary  and  independent  of  n(t).  The  statistics  of 
a(t)  must  be  chosen  so  that  y(t)  is  an  accurate  description  of  the  received 
noise.  We  will  further  make  the  reasonable  assumption  that  the  modulating 
process,  a(t),  is  slowly  varying  (compared  with  n(t))  so  that  the  spectrum  of 
a(t)  has  negligible  overlap  with  the  spectrum  of  n(t).  That  is,  n(t)  is  a 
bandpass  process,  while  a(t)  is  a  lowpass  process.  The  details  and 
arguments  as  to  why  this  is  a  reasonable  model  for  physical,  as  well  as 
mathematical  reasons,  are  given  by  Hall  (1966). 

An  appropriate  choice  for  a(t)  turns  out  to  be  a  process  with  a 
probability  density  function  given  by 


where  m  and  a  are  parameters  defining  a  two-sided  chi  distribution,  X2^n?»<r)- 
The  zero  mean  Gaussian  process,  n(t),  is  completely  described  by  its 
covariance  function,  R  (t  ),  and  its  probability  density  is  given  by 
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where  R  (0)  =  a  ,  the  variance  of  nit).  Since  n(t)  and  a(t)  are  assumed  to 
n  1 

be  statistically  independent. 


py<x,  =  f  FT  pa<!>  Pn«I> 


(A -4) 


X 

If  we  let  y  =  m^o^/cr  and  O  =  m  +  1  >  1,  then 


nf. 


©-  1 


pv(>0 


F^T 


<  X  <  CD  . 


(A-5) 


For  this  model  Hall  has  computed  the  distribution  of  the  noise 
envelope,  the  average  number  of  level  crossings  of  the  noise  envelope, 
and  the  distribution  of  level  crossings  (i.  e. ,  time  between  bursts)  and 
has  shown  that  by  proper  choice  of  the  parameters,  0  and  y,  the  proposed 
model  fits  quite  well  these  measured  statistics  for  LF  atmospheric  noise. 
Even  so,  there  are  some  problems.  First,  when  0  is  in  the  range 
2  <  0  £ 3,  y(t)  has  infinite  variance  (i.e. ,  infinite  energy)  and  therefore 
cannot  be  a  model  for  physical  noise,  even  though  it  is  found  to  fit  the 
measurements.  The  problem  arises  in  the  tails  of  the  distribution  of 
y(t)  being  such  that  y(t)  is  either  barely  convergent  or  barely  divergent, 
depending  on  the  parameter  0.  This  in  turn  arises  from  allowing  the 
modulating  process,  which  represents  the  statistics  of  the  noise  sources, 
to  have  an  infinite  range.  Note  that  y  and  0  can  be  related  to  the  two 
parameters  commonly  used  to  define  the  distribution  of  the  noise  envelope 
for  atmospheric  noise,  V  ,  and  the  rms  level,  F  (Spaulding  et  al. ,  1962). 
As  we  shall  see,  the  problem  of  infinite  variance  will  not  bother  us, 
because  we  will  simply  normalize  our  analysis  to  the  point  on  the  envelope 
distribution  that  corresponds  to  the  actual  rms  level  measured.  That  is, 
instead  of  normalizing  to  the  (sometimes  nonexistent)  rms  level  of  the 
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model,  we  will  use  a  level  that  corresponds  to  the  rms  level  of  the  actual 
measured  atmospheric  noise. 

If  we  are  interested  in  HF  atmospheric  or  man-made  noise,  we 
can  modify  the  above  model  as  follow'.:  since  HF  atmospheric  noise  and 
man-made  noise  generally  exhibit  a  much  smaller  dynamic  range  than 
LF  atmospherics,  the  modulation  process,  a(t),  can  be  allowed  to  have 
only  a  finite  range,  -  0<  a  <3.  While  this  introduces  another  parameter, 
3,  and  complicates  the  mathematics,  this  not  only  results  in  finite 
energy  for  y(t)  (for  all  choices  of  0),  but  should  make  it  possible  to 
model  man-made  noise  quite  closely  by  proper  specification  of  covariance 
functions  for  a(t)  and  n(t).  That  is,  we  will  use  for  the  distribution  of  the 
process  a(t). 


Pa<X>*  7"^Tl 

IX1 


(— 1 %  -  P5*  5  P  • 

Zv  x 


(A-6) 


where  k  is  chosen  so  that 


^  pa(x)dx  =  1 


Figure  1  shows  the  APD's  (envelope  distributions)  for  a  range  of 
values,  0,  and  a  value  of  y  that  corresponds  to  the  expected  APD's  in 
1000  Hz  bandwidth  in  the  northwest  United  States.  Figure  2  shows  the 
kind  of  modification  that  results  in  the  envelope  distribution  (0  =  3)  when 
various  values  of  P  are  used.  The  expressions  for  the  envelope  distribu¬ 
tions  have  been  derived  by  Hall  (1966)  and  are  given  in  appendix  B.  The 
APD's  shown  in  figures  1  and  2  were  calculated  by  computer. 


A.  3  Signal  Representation  and  Determination  of  the  Optimum  Receiver 

l  '  « 

In  this  section  we  will  derive  the  optimum  receiver  for  noise 

t 


and  signal  situations  pertinent  to  the  present  problem.  As  shown  by 


»**««»*' »*»* 


time  duration  of  the  signals  must  be  many  times  the  reciprocal  of  the 
receiver  RF  bandwidth.  In  fact.  Hall  has  shown  that  the  linear  matched 
filter  receiver  is  optimum  when  the  signal  duration  is  on  the  same  order 
as  the  receiver  reciprocal  bandwidth.  This  is  analogous  to  current 
noise -limiting  techniques,  in  which  large  bandwidths  (compared  with  signal 
duration)  are  required.  Typically,  we  may  use  a  signal  length  of 
0.01  sec  (100  bpsj/sec  with  an  RF  bandwidth  of  1000  Hz,  so  that  the  signal 
length  is  10  times  the  reciprocal  bandwidth. 

Let  us  consider  the  problem  of  deciding  which  signal  from 
among  a  set  of  two  signals,  s  (t)  and  s  (t),  is  represented  by  some 

X  Cd 

received  waveform,  x(t).  We  want  to  develop  a  decision  scheme  that 
will  minimize  the  average  "risk. "  The  "risk"  involved  in  making  an 
error  is  defined  as  the  probability  of  making  that  error  times  the  "cost" 
of  making  that  error.  The  problem  is  how  to  choose  between  two  hypotheses. 


H^  and  H^,  where 


^  :  x(t)  =  Sj(t)  +  a(t)  n(t) 

H2  :  x(t)  =  s2(t)  +  a(t)  n(t)  . 


(A -7) 


We  denote  the  "cost"  associated  with  H^  by  (the  cost  to  the 
observer  if  he  decides  H2  is  true  when,  in  fact,  is  true  minus  the  cost 
to  the  observer  if  he  decides  H^  is  true  when,  indeed,  it  is  true)  and  the 
cost  associated  with  H  by  C  .  If  the  probability  of  occurrence  of  s  is 

id  id  X 

q  ,  and  the  probability  of  occurrence  of  s  is  q  ,  we  have  the  threshold, 

A  id  id 

Q,  where 


\cx 

->2C2 


(A-8) 
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Let  L[x(t)  j  denote  the  likelihood  ratio,  given  by 

pUlt)^) 

=  p(x(t)|H2) 

The  minimum  average  risk  is  achieved  then  by  choosing  when 

L[x(t)]  >  Q  , 


(A-9) 


(A -10) 


and  choosing  H  otherwise. 
c • 


The  probability  of  making  the  error,  "decide  s  when  s  was  sent," 

X  4 


is  then 


P  =  prob  { L[x(t)]  >  Q| s  }  , 

C  Cm 


(A -11) 


so  that  the  probability  density  of  the  likelihood  ratio  is  required. 

In  the  completely  general  case  given  by  (A -7),  the  probabilities 
in  the  likelihood  ratio  are  extremely  difficult  (if  not  impossible)  to  compute. 
In  the  situations  of  interest  here,  in  which  the  receiver  bandwidth  is  sub¬ 
stantially  less  than  the  band  center  frequency,  the  hypotheses  (A -7)  can  be 
formulated  in  terms  of  the  complex  envelope  of  x(t)  (Helstron,  I960). 

Without  loss  of  generality,  we  can  lfet  Q  =  1  (equal  probabilities  of  occurrence 
and  equal  costs  associated  with  s  and  s  ),  and  we  will  further  let  s  (t)  =  0. 

X  u  b 

Hall  has  evaluated  the  likelihood  ratio  for  the  case  m  =  1.  We  will  follow 
Hall's  procedure  but  will  carry  the  calculations  through  for  arbitrary  m. 

Our  two  hypotheses  can  now  be  given  in  the  following  form: 


Hj  :  x(t)  =  p(t)  +  T](t) 


H2  :  x(t)  =  TXt) 


OStsT  , 


(A-12) 


Where  T  is  the  time  duration  of  the  signal,  p(t)  denotes  the  complex 
envelope  of  the  signal,  T)(t)  denotes  the  complex  envelope  of  the  additive 
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noise  process,  y(t),  and  x(t)  denotes  the  complex  envelope  of  the  received 
narrowband  waveform,  e.g.  , 

x(t)  =  x  \t)  cos  ij)  t  -  x  (t)  sin  u  t 
c  os  o 

=  R  (x(  , 

e 


X<t)  =  xc(t)  +  i  x#(t)  . 

Let  the  signal,  p(t),  be  represented  by  a  N  dimensional  vector,  p,  i^e. , 

jj>  =  {nj,  •  •  • »  wl*ere  the  p.  are  complex  numbers  and  are  samples 

from  the  waveform  ^(t).  We  will  now  assume  that  the  bandwidth  of  the 

signal  is  substantially  less  than  the  receiver  bandwidth,  2B.  The  signal 

may  then  be  represented  by  2BT  samples.  We  break  the  interval,  [o,  T], 

into  N  subintervals.  At.  and  n.  =  ^(t. ),  where  t.  denotes  some  time  in  the 

1  i  i  i 

interval.At..  The  waveform,  p(t),  can  be  represented  exactly  in  this  manner. 

Let  the  complex  noise  process,  TKt;,  be  also  represented  by  N 
samples,  30  that  T).  =  T](t.  )(At.)^.  The  "H.  are  now  complex  random  variables, 
and  the  process,  T|(t),  is  represented  by  the  random  vector  ][=  {T^,  T)  ,  . . . ,  TJ^j. 
(For  the  problems  involved  in  representing  a  random  waveform  by  a  finite 
dimensional  random  vector  and  the  connection  between  the  processing  of 
random  vectors  and  the  processing  of  continuous  random  waveforms,  see 
Hannock  and  Wintz  (1966,  app.  B). 

Because  of  the  above  constraint  on  the  signal,  we  may  assume  that 
the  various  T)  are  independent  (which  corresponds  to  the  atmospheric 
noise  being  white  prior  to  detection). 

Our  two  hypotheses  are  now  given  in  vector  form: 


H1  :  X  =  JJ  +  ][ 

H2  :X  =  1 


and  we  must  now  compute  p(x  { H  )  and  p(x  |  H  )  . 

X  L> 


(A-13) 
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pix  i  Hj)  =  (X  "  £)  . 


(A -14) 


so  that  (T))  is  required.  Recalling  that  the  noise  process,  y(t),  is 
given  by  a(t)n(t),  we  can  write 


\  =  a.  v.  , 
i  11 


(A-15) 


where  a  =  a(t.)  and  v.  =  v(t.)  (At.)  ,  in  which  v(t)  is  the  complex  envelope 
1  1  i  i  i 

of  n(t).  Since  n(t)  is  a  Gaussian  process,  v(t)  is  a  complex  Gaussian 
process.  Thus  the  probability  density  function,  p^  (v),  of  the  vector,  v, 
can  be  written  in  the  form 


pv  (->  =  "fTTr  ]  . 

"1*1 


(A  -16) 


where  v  denotes  the  conjugate  of  v,  and  v_  ,j,  denotes  the  transpose  of 
the  vector,  v  ;  $  is  the  N  by  N  covariance  matrix  of  v(t),  i.  e. , 


$=  {E(v.  v.)},  1  s  i»  j  s  N,  where 

E(v*  v.)  =  E[  Mt.)  v(y]  (At.Ay  * 

where  E  denotes  the  expectation  operator. 

Our  assumption  of  impedance  states  that 

E  l  v*(t4)  v(V)]  =  <Tj26„,  i,  j  =  1,  2.  •  •  • »  N 


(A -17) 


(A -18) 


where  6..  is  Kronecker’s  delta,  i.  e. ,  6. .  =  0  for  i  /  j,  and  6. .  =  1,  i  =  j. 
ij  ij  ij 


Letting  AL  »  At  for  all  i,  we  obtain  then 

Tcbr  *xp[\~z7~  -T^J  ■ 


IT  {‘"j  At) 


ri  At 


(A -19) 


To  compute  p^tjl)  we  will  make  the  transformation  b  =  i/a,  where  now, 
since  the  distribution  of  a  is  given  by  (A -2),  the  distribution  for  b  is 


m  Xn/Z 

Pb(b)=  ~m“'m  lb|m“1*xp[-  -^b2]  . 

•  n-g-)  L  2^  J 
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(A-20) 


and  then 


( A-21) 


PJ1-  =  Pb*-  d-  * 

where  b  is  the  vector  {b,,  b  ,  ...b  },  and  b.  =  l/a..  Our  assumption  of 
—  1  2  n  ii 

s  independent,  i.  e., 

E[b{t.)b(t.)]  =  <r26..  . 


independence  makes  the  b/s  independent,  i.e.. 


(A-22) 


Therefore  the  distribution  of  the  vector,  b,  based  on  (A -20) 

„  m/2 

m  ' 


p  ,y  .  -2_ 


iUi  lbi!m‘‘exi>[-71  tTk] 

2<r 


(A -23) 


Now,  since  T],  =  v./b.,  I]  and  v  complex,  the  Jacobian  of  the  transformation 

1,11 

v.  -*  71.  ib'  b.".  This  means,  then, 

4  11 


pa(a|-)  =(11I1  bi2)pv(,,i *  biV  1  *  >• 2 . 


N)  . 


Therefore  from  (A-19), 

p^aiy^Tijh2) 

Substituting  into  (A-21) 


p  (n)  =  K 
R~ 


N 

1  F 

1  } 

*N(a 

2a  ,N  expL 

1  At)  L 

ffj2  At  i  =  1 

N 

IT 

, ,  ,  m  +1 
b  .  eXp 

[-l(~ 

i  =  l 

*  ! 

L  i  =  l  V  2<r 

(A-24) 


»  (A-25) 


where 


K  = 


<T> 


m/2 


~rrr 


m  <X)  N  2  ,N 
(T  F  IT  (ffj  At) 


(A-26) 


Evaluating  (A-25),  we  obtain 

N 


m  +  2 


P  (ti)  -•  KTl  r  fllLli')  ( JjliL  +  _£EL  ^ 

a-  i  =  i  A  2Af  + jy 


(A-27) 


“•j  At  2cr 
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m  +  2 


?.  n 


jj,( I"/*"**-)  • 


From  (A. -14),  our  likelihood  ratio,  L[x],  becomes 

N 


m  +  2 
ZA*  “Z 

mcr,  At 


Mx)  = 


r  ■»  mtrt  At  *i 

i?l[lxil  +-^T—  3 


m  +  2 


Nr  7  mcr^At  ] 


Our  optirm,:ti  decision  rule  (A -4)  becomes  "choose  Hj,  if 


N 

n 


i  =  lU  * 


e. 

mcr  c*  t 
!  X- !  *  +  — K — 


2<r 


4:[lVM 


2  ii 

2  mcr  At  ■} 


This  is  equivalent  to  the  rule 
N  2 


i  =  l  i  =  1 


(A-28) 


(A-Z9) 


(A-30> 


(A-31) 


and  in  terms  of  the  time  functions  available  to  the  receiver,  the  rule 
(A-31)  becomes 


f  in  [|X(t)|2  + 
0 


|X(t)  -  HL(t)|  2  + 


dt,  (A-32) 


There  are  many  physical  realizations  of  (A-32)  (see,  for  example.  Hall, 

1966).  The  receiver  consists  of  two  branches,  one  to  compute  each  side  of 

the  inequality  (A-32).  Each  branch  could  contain  an  envelope  detector, 

followed  by  a,  square  law  device  (which  gives  |x(t)l  ,  for  example), 

2  2 

followed  by  a  log  amplifier  biased  by  mcr^  /2a  (a  bias  proportional  to  the 


-89- 


actual  expected  noise  power),  followed  by  an  integrator.  Threshold  logic 
circuitry  would  then  compare  the  thus  computed  quantities  and  announce 
a  decision.  This  receiver  can,  therefore,  be  easily  constructed.  Any 
digital  receiver  is  simply  a  box  for  computing  some  quantity  from  the 
received  waveform  and  for  comparing  this  quantity  with  a  threshold. 

Figure  3  (from  Hall ,  1966  )  shows  a  block  diagram  of  one  receiver  structure 
which  implements  the  decision  rule.  It  now  remains  to  analyze  the  per¬ 
formance  (compute  the  error  rate)  of  our  receiver. 

A. 4  Calculation  of  Probability  of  Error 
The  probability  of  error  expression  (A -11)  is 


p  =  prob  {L(x)  >  1  I  H  is  true}  , 
or,  from  (A -29),  since  being  true  implies  |x(t)|  =  |Tl(t)|  , 


(A-33) 


p^  =  prob 


r  N 

{ n 
^  i  =  1 


2 

m«r  At 

l”il  *1^— 

2 

m<r  At 

IVM  ~ 


(A-34) 


Putting  (A-34)  in  a  more  manageable  form,  we  have 


=  prob  ■^^fnZj>o|‘  , 


(A-35) 


where 


z.  = 

i 


2 

me.  At 
2  1 
IM  + 

__  ^ 

mcr^  At 

'vMz  +  — 


(A-36) 


Now,  it  is  not  possible  to  compute  o  directly,  since  we  cannot 

N  ® 


obtain  the  probability  density  for 


1  lnzi 


,  and  we  must  be  satisfied  with 
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good  estimates  for  p  .  Since  p  is  given  in  terms  of  a  sum  of  random 

0  0 

variables,  the  characteristic  functions  for  the  random  variables,  Inz., 
should  be  useful.  We  will  use  characteristic  functions  to  obtain  the 
appropriate  Chernoff  bound  on  p^  (Van  Trees,  1968). 

Let  4>(s)  denote  the  characteristic  function  for  our  likelihood 
ratio  L,  under  hypothesis  H^,  that  is, 


4>(s)  =  f 


eSL  P(L|H  )  dL  . 


Now  define 


(A-37) 


p(s)  =  I  n<j>(s)  . 


We  then  have  the  Chernoff  upper  bound 


Pe  1 


<eP(»)  -  ®P‘S>,  ,>0  . 


(A-38) 


(A  -39) 


and  the  minimum  upper  bound  is  obtained  by  evaluating  (A -39)  for  the 
value  of  s,  sq,  which  makes  p(s)  =  0.  We  can  now  apply  the  central 
limit  theorem  to  the  above  to  obtain  the  approximation  (Van  Trees,  1968) 


Pe  =i{exp(p(s)  -  sp(s)  +  —  p  ( s )] }  erfc  [s  n/  f>*(s)  , 

where  erfc  denotes  the  complimentary  error  function.  For  values  of 
s V  p'(  s )  >  3,  we  can  approximate  the  erfc  term  to  obtain 


(A -40) 


exp  [  p(s)  -  sp(s)] 


Ire  s2p(s) 


(A- 41) 


It  turns  out  that  for  the  range  of  pg  of  interest  to  us  (10  and  less)  the 

factor  1/'^2tts^(T)  no*  aPPrec^a^^y  change  the  results,  and  besides, 

•  • 

as  we  shall  see,  it  is  extremely  difficult  to  evaluate  p(s  ). 

o 

In  our  case,  then,  if  g^(  s)  denotes  the  characteristic  functions 
for  the  independent  random  variables  I nz^, 


p  ( s )  =  f  n  T1  g.(s)  , 

i=  1  1 


(A-42) 
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so  that  the  approximation  we  will  use  is  given  by 

It  follows  that  (Hall,  1966) 


(A-43) 


*.(•)  =  E[exp  (sfn  z.)] 


n 

r«f> 


f J  -s.  j 


dn. 


(2s  -  m  -  ?)/2 


is 


where 


f . 


2<r 


2  ,  and  denote  the  real  and  imaginary  parts 


2a 

mo-j  At 


of  the  complex  number  T)  .  Equation  (A-43)  states  that  the  performance  is 
dependent  on  the  actual  noise  power  (which  is  proportional  to  §^)  and  on 
the  actual  signal  shape  and  energy  (i.e.,  depends  on  the  2TB  samples, 
y .,  of  the  signal  complex  envelope).  To  simplify  things,  let  the 
imaginary  part  of  y.  =  0,  for  all  i.  ,An  example  of  such  a  signal  would  be 


•£t)  =  r(t)  co*  V  • 


The  samples,  y.,  are  now  real  numbers,  and  our  integral  (A-44)  reduces  to 


2s  -  m  -  2 


<*)  = 


r( 


m  +  2 


) 


fdx  r  dy  & 

*  ^  [{x -y  )2+ y2  + 


r(~)  w 


(A-45) 


ft 


I 


I  4 


i 


As  one  might  expect,  and  as  Hall  (1966)  has  shown,  the  best  signal  design 
is  a  signal  whose  energy  is  uniformly  distributed  in  the  signaling  interval. 
This  means  that  the  y.  can  be  considered  equal,  resulting  in  needing 
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i 


to  evaluate  only  one  g.(s).  Our  estimate  now  becomes 


N 


Pei1T  gi<s0>=[g(s0>] 

i  =  l 


2BT 


(A-46) 


In  terms  of  the  total  energy  per  received  bit,  and  the  noise  power 
density,  (A-45)  becomes 

2s  -  m  -2 

r(iH-±A) «,  .  2  2  2 

g(s)  -  -  fdx  Pdy  — — *  y  *  1 _ 

*r(4r-)  J  -2  2  7  9 

^  **  (*  +y  -2ax  +  a  +1) 


(A -47) 


where  the  parameter,  a,  is  given  by 


•i  * 


l 

i 


•Ttbn" 


(A -48) 


where  E  denotes  the  signal  energy  (joules)  and  N  denotes  the  noise  power 
density  ( watts/hertz  ).  The  integral  (A-41)  must  be  evaluated  for  8^80» 
where,  from  (A -42),  sq  is  the  solution  to 


f.J* 


(x2  +  y2+l) 


2s  -  m  -  2)/2 


2  2  2 
(x  +y  -  2ax  +  a  +1) 


<*44— rr) 


to  to 

M 


dy 


(x2  +y2  +1) 


x"  +  y”  -  2 ax  +  a“  +  2  '  (A-49) 

(2s  -  m  -  Z)/Z 


-•  -w  (x2  +  y2  -  2ax  +  a2+l) 


For  each  a  of  interest,  (A-49)  must  be  solved  for  s  and  then  evaluated 

o 

for  this  value  of  s.  Quite  fortunately  it  can  be  shown  that  (A-49)  is 
uniquely  solved,  for  all  a,  by 


8 

O 


m  4  2 
4 


(A-50) 
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The  integral  (A-47)  must  now  be  evaluated  for  this  sQ»  and  computer  tech¬ 
niques  have  been  developed  and  programs  written  to  do  this. 

The  case,  0=3  (m  =  2),  is  a  special  case.  For  m  =  2,  (A-49) 
is  uniquely  solved  for  all  a  by  s  =  1,  and  for  s  =  1,  the  integral  (A-47) 
can  be  evaluated  in  closed  form  by  first  performing  the  x  integration  and 
then  the  y  integration.  The  result  is 


g(l>=  — -  tn(-£*±±±)  . 

a^a2  +4  W  «2  +  4  -  c 


(A -51) 


From  (A -46),  we  see  that  we  can  make  the  probability  of  error 
arbitrarily  small  by  making  2BT  large  enough.  Of  course,  it  is  well 
known  that  by  signaling  slowly  enough  (i.  e. ,  making  T  large),  the 
probability  of  error  can  be  made  as  small  as  desired. 

We  will  now  consider  the  following  example:  summer  nighttime, 
Washington  state  area,  <Dq  =  500  kHz,  1000  Hz  bandwidth,  and 
100(2TB  =  10)  and  50  (2TB  =  20)  bps.  From  CCIR  (1964)  we  have  F 

am 

F  (500  kHz)  =  133  dB,  and  V,  (1  kHz  bandwidth)  =  14  dB,  i.e. , 
am  dm  ^ 

dynamic  range  of  the  envelope  distribution  for  probabilities  from  10 

to  .99  is  about  88  dB.  This  corresponds  to  0  =  3.  F  =  133  dB  gives 

8  s.m 

us  that  the  noise  power  density  is  7.95  x  10  W/  Hz.  This  number  is 
required  to  determine  the  value  of  yto  use  in  specifying  the  envelope 
distribution  and  the  amount  of  bias  required  on  the  log  amplifiers  of  our 
optimum  receiver. 

Figure  4  shows  the  performance  of  the  optimum  receiver.  It  also 
shows  the  performance  of  a  binary  matched  filter,  differentially  coherent 
phase  shift  keying  receiver  for  both  Gaussian  noise  and  atmospheric 
noise  (0  =  3). 
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If  we  now  tried  to  improve  our  estimate  by  using  (A-40)  or 

•  • 

(A -41),  we  would  need  to  evaluate  p  (s)  for  s  =  1,  m  =  2.  This  is 

2TB ^dx  f  dy  [in  JL+Y**  1 

p(l)  -  ~'a>  (x  -  a)  ty  +1  L  (x  -  a)  -f  y  +  1  J 

oo  at  2  2 

U  fdy  ty, 

-«  -S  (x  -  a)  +y  +1 


(A-52) 


However,  instead  of  trying  to  evaluate  the  above,  we  note  from 
figure  4  that,  if  the  factor  IA/^tTs^p^s)  were  on  or<^er  °*  0.1,  this 
would  change  our  estimate  by  less  than  1  dB,  and  even  if  1/V2tt  s^'p (s) 
were  on  the  order  of  0.  01,  a  change  of  less  than  2  dB  would  result. 

That  is,  for  a  given  probability  of  error,  the  required  SNR  would 
be  not  more  than  2  dB  less.  This  indicates  that,  while  the  curve  shown 
is  an  upper  bound  on  the  true  probability  of  error  characteristic,  it  is 
a  very  tight  upper  bound.  Estimates  of  (A-52)  indicate  that  the  factor 
1/v2n  s2  *P(s)  *s  somewhere  between  0.1  and  0.25. 

Figure  5  shows  the  performance  of  the  optimum  receiver  for  a 
range  of  values  of©  for  2  TB  =  10,'  and  figure  6  shows  the  same  for 
2TB  =  20. 

If  we  try  to  carry  out  the  above  calculations  with  a  modified 
Hall  model  (to  fit  HF  atmospherics  or  man-made  noise),  we  find  that, 
besides  being  much  more  difficult,  the  optimum  receiver  is  different 
for  each  value  of  3  we  might  use.  Studies  like  this  with  a  modified 
model  remain  to  be  done. 

For  completeness,  we  have  included  figures  7,  8,  and  9,  which 
summarize  the  performance  of  typical  linear  receivers  in  noise  of 
atmospheric  type.  Figure  7  shows  the  relationship  between  the  proba¬ 
bility  of  error  and  the  signal -to -noise  ratio,  parametric  in  V  ,,  for 

Cl 


binary  differentially  coherent  phase  shift  keying  (DCPSK).  These  curves 
were  obtained  by  the  method  described  by  Halton  and  Spaulding  (1966) 
and  also  are  applicable  to  noncoherent  frequency  shift  keying  (NCFSK), 
if  we  add  3  dB  to  the  abscissa  value. 

Figure  8  shows  the  probability  of  element  error  for  a  four -phase 
DCPSK  system  as  a  function  of  signal -to -noise  ratio,  parametric  inV^. 
As  in  figure  7,  the  signal  is  assumed  to  be  steady.  The  binary  error  rate 
is  approximately  equal  to  one -half  the  element  error  rate  at  large 
signal -to -noise  ratios. 

Figure  9  shows  probability  of  error  characteristics  for  both 
two-  and  four -phase  DCPSK  systems  when  the  signal  is  Rayleigh  fading, 
as  described  by  Halton  and  Sjpaulding  (1966). 
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Envelope  Level,  dB>  y 
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Percent  of  Time  Ordinate  is  Exceeded 

Figure  A.  2  Amplitude  probability  distributions  of  the  noise  envelope  of  the 
modified  Hall  model  for  a  range  of  values  of  the  parameter  3. 
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APPENDIX  B 


Distribution  of  the  Envelope  and  Phase  for  the  Hall  Model 


Our  narrowband  noise  process  is  of  the  form 

y(t)  =  a(t)n(t)  . 

Since  this  is  a  narrowbard  process,  it  can  be  expressed  as 


(B  -1) 


[«t  +  4>(t)  , 


y(t)  =  V(t)  cos  jw^t  +4>(t) 

where  V(t)  is  the  noise  envelope  process, and  <j>{t)  is  the  phase  process. 
The  joint  probability  density  of  V  and  is  given  by 


(B-2) 


P„  ,(V,  <J))  =  ’/p  ~  (Vcos<j>,  Vsin<|>)  , 

V,4>  y»  y 


(B-3) 


where  y  denotes  the  Hilbert  transform  (quadrature  component)  of  y. 

That  is,  the  envelope  V  is  given  by 

l 

V=(y2  +  y2)2.  (B-4) 

Now  since  the  spectrum  of  a(t)  has  negligible  overlap  with  the  spectrum 
of  n(t),  we  have 

y(t)  =  a(t)  n(t)  , 

(B-d) 

and  since  n(t)  is  a  Gaussian  process,  n(t)  is  a  Gaussian  process,  indepen¬ 
dent  of  n(t),  and  identically  distributed  as  n(t). 

Thus, 


py.  y(y'  Y>  ‘  pan,  r> 


P  -{x, 
a.n.n  x  x 


-00 


Pa(x>  Pn^>  pi$ 


(B-6) 


-106- 


and 


P„<x)  «  P~M  «  (Zv,*)  , 

2c. 


and 


Pa(x)  « - £ 

a  ,  .  m 


+  1  exp(Ti  “T/  '  <  p. 


2fr  X 


where  k  is  chosen  to  satisfy 


I 

-P 


p  (x)  dx  =  1  . 

a 


So,  we  have.  Substituting  (B-7)  and  (B-8)  into  (B-6), 


py  ~(y.  Y)  = 


Substituting  into  (B-3)  gives 


vk  r  dx 


k  i 

P 

"  dx  / 

£  m  y2  y2  "N 

7  ' 

ir«r.£ 

1  f 

)  xm+  3  p 

1 

^2cr^x2  2(r2x2  2,r2x2^ 

■ 


Py  =  ~z~  \  - - -  exp  1 

*  *  t  2  }  m  +  3  *  \  ^  i 

TT  (T  v  y  \  c  -4  4‘  2  2 

1  0  Za  x  2^  x  2tr  x 


f  _  -m  .  V2cosZ4>  V2sin2< 
V  -2  2  .'22  ~  7  7 


P 

_  _Vk  n  dx 


“p(\.2  z 

1  0 


m 


2o"”x“  2cTj2xZ 


')• 


1 


Let  z  -  ~  ,  then 
x 


V,4> 


&) 

Vk  r  m/2  f 

2  J  Z  exP  1 

*1  l/f* 

m 


2«Z  2„,2 


')] 


dz 


So,  we  see  that  <j>  is  uniformly  distributed,  i.e.. 


P<|>(<i>)  *  27  *  0  <  *  <  2^ 


(B-7) 

(B-8) 

•  (B-9) 

(B-iO) 

■) 

(B-n) 

(B-12) 
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"A  w 


and 


p  {v)s™ 

pv'  '  2 

^  1/P2 


dz 


(B  -13) 


The  integral,  Py(V)  is  easily  evaluated  in  closed  form  for  m  even,  and 
is  easily  evaluated  numerically  (by  Gauss  -Laquerre  quadratures)  for 
any  m. 

For  p=  w  ,  we  have 


pv<v>  = 


m 

Y  V 


m  +  2 


(B-14) 


(v2*v2) 


and  for  the  probability  of  the  envelope  exceeding  the  level  V^,  that  is, 
the  APD, 


P(V  >  Vq)  =  f  py(x)  dx  = 


Z  m 

[(?>■] 


72  * 


( B  —15) 
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APPENDIX  C 


Average  Rate  of  Envelope  Level  Crossings 


Letting  N(Vq)  denote  the  average  number  of  crossings  per  second 
of  the  level  Vq  by  the  noise  envelope,  V,  it  is  well  known  (Middleton,  1961, 


p.  427)  that  for  a  stationary  random  process,  V(t), 

N(Vo)  =  J|V|[pViV(v,V)]dV  . 


(C-1) 


where  V(t)  =  dV(t)/dt. 

As  shown  in  appendix  A,  the  assumption  of  negligible  spectral  over¬ 
lap  of  a(t)  and  n(t)  allows  the  envelope  of  the  noise  process  to  be  given 


V(t)  =  |a(t)  |E(t)  , 


(C-2) 


where  E(t)  is  the  envelope  of  the  narrow -band  Gaussian  process  n(t). 


Thus, 


,  Ut)1  +  E(t) 


(C-3) 


If  we  are  dealing  with  bandwidth s  that  are  a  small  fraction  of  our 
center  frequency,  uj^,  it  is  reasonable  to  assume  that  the  modulating 
process,  a(t),  is  varying  sufficiently  slowly,  so  that 


(C-4) 


which  allows  us  to  represent  V(t),  where 


V(t)  =  |a(t)|E(t). 


(C-5) 


PV.V(V'V>  =  P,a|E,,a,E(V^ 


(C-6) 


-109  - 


rtWWBMffWf' 


i 


where,  from  appendix  B, 

/  *  2k 


—  exf(-7f t)-  0 


<  x  <  (3  . 


(C-7) 


2cr  x 


The  problem  now  reduces  to  finding  p^  •  (E,E).  Since  E(t)  is  the 
envelope  of  a  narrow -band  Gaussian  process,  it  can  be  shown  that,  in 
general  (Middleton,  1961,  p.  420), 


E  /'^33E 

PF  p(E»E)=  - 1  exP\. - 

(2-  )*  V  2|pJ 


^  r'EWj3^2i4>  i  , 

>xpl — 53m - J  '  |C-f 


where  |n|  « (bbj  -  1^  )  .  =  »>2|^p  .  I*jj* 


and  the  b^  are  the  frequency  moments 


d  =  fs  I 
n  J  n 


(f)(w  -  wq)  df  , 


(C-9) 


where  S  (f)  denotes  the  spectral  density  of  the  process  n(t).  Since  S  (f) 
n  n 

is  the  Fourier  transform  of  the  correlation  function  for  n(t),  R  (x  ), 

n 

these  moments  can  also  be  computed  from  R^t  )• 

If  the  spectrum  is  symmetrical  about  then  (C-8)  reduces  to 


PE..E(E‘E)=b  sfa bl 

O  Z 


exp[-fe-  -  sr]  •  <CJ0) 

t-  n  7  J 


Then,  with  z  =  l/x  , 


pv,  v(v*  v)  = 


V/2vb2  j^p2 


m  +  1 
2 


[-(-*. !♦  —  +  -£)] 

L  '2<r  2bc  2bz/J 


dz  .  (C-ll) 


The  frequency  moments,  b^,  are  determined  by  the  bandpass  shape  of 
the  receiver,  and  once  known,  then  N(VQ)  can  be  ev^^ted. 
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Now,  if  R  (t  )  has  the  typical  form 
n 

R  (t  )  =  R  (t  )  cos  to  t  , 
n  c  o 


and  vve  define  the  rms  bandwidth.  B_, 


R  (o) 
c 


(C-13) 


where  S  (f)  is  the  Fourier  transform  of  R  (t  ),  Hall(1966)  has  shown  that, 
c  c 


when  (3  =  co  , 


N(V  =  -  - - <- - ---  , 

°  /e-iW.  2.  zW2 


°  r(^)eVV+v2) 

l/2 

where,  as  before,  y  =  m  '  <^/<r  and  ©  =  m  +  1  >  1. 


( C  -14 ) 
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